Addressing the batch effect issue for LC/MS metabolomics data in data preprocessing.

With the growth of metabolomics research, more and more studies are conducted on large numbers of samples. Due to technical limitations of the Liquid Chromatography-Mass Spectrometry (LC/MS) platform, samples often need to be processed in multiple batches. Across different batches, we often observe differences in data characteristics. In this work, we specifically focus on data generated in multiple batches on the same LC/MS machinery. Traditional preprocessing methods treat all samples as a single group. Such practice can result in errors in the alignment of peaks, which cannot be corrected by post hoc application of batch effect correction methods. In this work, we developed a new approach that address the batch effect issue in the preprocessing stage, resulting in better peak detection, alignment and quantification. It can be combined with down-stream batch effect correction methods to further correct for between-batch intensity differences. The method is implemented in the existing workflow of the apLCMS platform. Analyzing data with multiple batches, both generated from standardized quality control (QC) plasma samples and from real biological studies, the new method resulted in feature tables with better consistency, as well as better down-stream analysis results. The method can be a useful addition to the tools available for large studies involving multiple batches. The method is available as part of the apLCMS package. Download link and instructions are at https://mypage.cuhk.edu.cn/academics/yutianwei/apLCMS/

Particulate metal exposures induce plasma metabolome changes in a commuter panel study

Introduction Advances in liquid chromatography-mass spectrometry (LC-MS) have enabled high-resolution metabolomics (HRM) to emerge as a sensitive tool for measuring environmental exposures and corresponding biological response. Using measurements collected as part of a large, panel-based study of car commuters, the current analysis examines in-vehicle air pollution concentrations, targeted inflammatory biomarker levels, and metabolomic profiles to trace potential metabolic perturbations associated with on-road traffic exposures. Methods A 60-person panel of adults participated in a crossover study, where each participant conducted a highway commute and randomized to either a side-street commute or clinic exposure session. In addition to in-vehicle exposure characterizations, participants contributed pre- and post-exposure dried blood spots for 2-hr changes in targeted proinflammatory and vascular injury biomarkers and 10-hr changes in the plasma metabolome. Samples were analyzed on a Thermo QExactive MS system in positive and negative electrospray ionization (ESI) mode. Data were processed and analyzed in R using apLCMS, xMSanalyzer, and limma. Features associated with environmental exposures or biological endpoints were identified with a linear mixed effects model and annotated through human metabolic pathway analysis in mummichog. Results HRM detected 10-hr perturbations in 110 features associated with in-vehicle, particulate metal exposures (Al, Pb, and Fe) which reflect changes in arachidonic acid, leukotriene, and tryptophan metabolism. Two-hour changes in proinflammatory biomarkers hs-CRP, IL-6, IL-8, and IL-1β were also associated with 10-hr changes in the plasma metabolome, suggesting diverse amino acid, leukotriene, and antioxidant metabolism effects. A putatively identified metabolite, 20-OH-LTB4, decreased after in-vehicle exposure to particulate metals, suggesting a subclinical immune response. Conclusions Acute exposures to traffic-related air pollutants are associated with broad inflammatory response, including several traditional markers of inflammation

Plasma Metabolomics in Human Pulmonary Tuberculosis Disease: A Pilot Study

We aimed to characterize metabolites during tuberculosis (TB) disease and identify new pathophysiologic pathways involved in infection as well as biomarkers of TB onset, progression and resolution. Such data may inform development of new anti-tuberculosis drugs. Plasma samples from adults with newly diagnosed pulmonary TB disease and their matched, asymptomatic, sputum culture-negative household contacts were analyzed using liquid chromatography high-resolution mass spectrometry (LC-MS) to identify metabolites. Statistical and bioinformatics methods were used to select accurate mass/charge (m/z) ions that were significantly different between the two groups at a false discovery rate (FDR) of q<0.05. Two-way hierarchical cluster analysis (HCA) was used to identify clusters of ions contributing to separation of cases and controls, and metabolomics databases were used to match these ions to known metabolites. Identity of specific D-series resolvins, glutamate and Mycobacterium tuberculosis (Mtb)-derived trehalose-6-mycolate was confirmed using LC-MS/MS analysis. Over 23,000 metabolites were detected in untargeted metabolomic analysis and 61 metabolites were significantly different between the two groups. HCA revealed 8 metabolite clusters containing metabolites largely upregulated in patients with TB disease, including anti-TB drugs, glutamate, choline derivatives, Mycobacterium tuberculosis-derived cell wall glycolipids (trehalose-6-mycolate and phosphatidylinositol) and pro-resolving lipid mediators of inflammation, known to stimulate resolution, efferocytosis and microbial killing. The resolvins were confirmed to be RvD1, aspirin-triggered RvD1, and RvD2. This study shows that high-resolution metabolomic analysis can differentiate patients with active TB disease from their asymptomatic household contacts. Specific metabolites upregulated in the plasma of patients with active TB disease, including Mtb-derived glycolipids and resolvins, have potential as biomarkers and may reveal pathways involved in TB disease pathogenesis and resolution

Novel Metabolic Markers for the Risk of Diabetes Development in American Indians

OBJECTIVETo identify novel metabolic markers for diabetes development in American Indians.RESEARCH DESIGN AND METHODSUsing an untargeted high-resolution liquid chromatography–mass spectrometry, we conducted metabolomics analysis of study participants who developed incident diabetes (n = 133) and those who did not (n = 298) from 2,117 normoglycemic American Indians followed for an average of 5.5 years in the Strong Heart Family Study. Relative abundances of metabolites were quantified in baseline fasting plasma of all 431 participants. Prospective association of each metabolite with risk of developing type 2 diabetes (T2D) was examined using logistic regression adjusting for established diabetes risk factors.RESULTSSeven metabolites (five known and two unknown) significantly predict the risk of T2D. Notably, one metabolite matching 2-hydroxybiphenyl was significantly associated with an increased risk of diabetes, whereas four metabolites matching PC (22:6/20:4), (3S)-7-hydroxy-2′,3′,4′,5′,8-pentamethoxyisoflavan, or tetrapeptides were significantly associated with decreased risk of diabetes. A multimarker score comprising all seven metabolites significantly improved risk prediction beyond established diabetes risk factors including BMI, fasting glucose, and insulin resistance.CONCLUSIONSThe findings suggest that these newly detected metabolites may represent novel prognostic markers of T2D in American Indians, a group suffering from a disproportionately high rate of T2D

Metabolic Profiles of Obesity in American Indians: The Strong Heart Family Study

The authors would like to thank the Strong Heart Study participants, Indian Health Service facilities, and participating tribal communities for their extraordinary cooperation and involvement, which has contributed to the success of the Strong Heart Study. The views expressed in this article are those of the authors and do not necessarily reflect those of the Indian Health Service.Obesity is a typical metabolic disorder resulting from the imbalance between energy intake and expenditure. American Indians suffer disproportionately high rates of obesity and diabetes. The goal of this study is to identify metabolic profiles of obesity in 431 normoglycemic American Indians participating in the Strong Heart Family Study. Using an untargeted liquid chromatography–mass spectrometry, we detected 1,364 distinct m/z features matched to known compounds in the current metabolomics databases. We conducted multivariate analysis to identify metabolic profiles for obesity, adjusting for standard obesity indicators. After adjusting for covariates and multiple testing, five metabolites were associated with body mass index and seven were associated with waist circumference. Of them, three were associated with both. Majority of the obesity-related metabolites belongs to lipids, e.g., fatty amides, sphingolipids, prenol lipids, and steroid derivatives. Other identified metabolites are amino acids or peptides. Of the nine identified metabolites, five metabolites (oleoylethanolamide, mannosyl-diinositol-phosphorylceramide, pristanic acid, glutamate, and kynurenine) have been previously implicated in obesity or its related pathways. Future studies are warranted to replicate these findings in larger populations or other ethnic groups.Yeshttp://www.plosone.org/static/editorial#pee

Transcriptome–metabolome wide association study (TMWAS) of maneb and paraquat neurotoxicity reveals network level interactions in toxicologic mechanism

A combination of the herbicide paraquat (PQ) and fungicide maneb (MB) has been linked to Parkinson's disease. Previous studies show that this involves an additive toxicity with at least two different mechanisms. However, detailed understanding of mixtures is often difficult to elucidate because of the multiple ways by which toxic agents can interact. In the present study, we used a combination of transcriptomics and metabolomics to investigate mechanisms of toxicity of PQ and MB in a neuroblastoma cell line. Conditions were studied with concentrations of PQ and MB that each individually caused 20% cell death and together caused 50% cell death. Transcriptomic and metabolomic samples were collected at time points prior to significant cell death. Statistical and bioinformatic methods were applied to the resulting 30,869 transcripts and 1358 metabolites. Results showed that MB significantly changed more transcripts and metabolites than PQ, and combined PQ + MB impacted more than MB alone. Transcriptome–metabolome-wide association study (TMWAS) showed that significantly changed transcripts and metabolites mapped to two network substructures, one associating with significant effects of MB and the other included features significantly associated with PQ + MB. The latter contained 4 clusters of genes and associated metabolites, with one containing genes for two cation transporters and a cation transporter regulatory protein also recognized as a pro-apoptotic protein. Other clusters included stress response genes and transporters linked to cytoprotective mechanisms. MB also had a significant network structure linked to cell proliferation. Together, the results show that the toxicologic mechanism of the combined neurotoxicity of PQ and MB involves network level interactions and that TMWAS provides an effective approach to investigate such complex mechanisms

Antagonistic Interactions in Mitochondria ROS Signaling Responses to Manganese

Antagonistic interaction refers to opposing beneficial and adverse signaling by a single agent. Understanding opposing signaling is important because pathologic outcomes can result from adverse causative agents or the failure of beneficial mechanisms. To test for opposing responses at a systems level, we used a transcriptome–metabolome-wide association study (TMWAS) with the rationale that metabolite changes provide a phenotypic readout of gene expression, and gene expression provides a phenotypic readout of signaling metabolites. We incorporated measures of mitochondrial oxidative stress (mtOx) and oxygen consumption rate (mtOCR) with TMWAS of cells with varied manganese (Mn) concentration and found that adverse neuroinflammatory signaling and fatty acid metabolism were connected to mtOx, while beneficial ion transport and neurotransmitter metabolism were connected to mtOCR. Each community contained opposing transcriptome–metabolome interactions, which were linked to biologic functions. The results show that antagonistic interaction is a generalized cell systems response to mitochondrial ROS signaling

Addressing the batch effect issue for LC/MS metabolomics data in data preprocessing

Abstract With the growth of metabolomics research, more and more studies are conducted on large numbers of samples. Due to technical limitations of the Liquid Chromatography–Mass Spectrometry (LC/MS) platform, samples often need to be processed in multiple batches. Across different batches, we often observe differences in data characteristics. In this work, we specifically focus on data generated in multiple batches on the same LC/MS machinery. Traditional preprocessing methods treat all samples as a single group. Such practice can result in errors in the alignment of peaks, which cannot be corrected by post hoc application of batch effect correction methods. In this work, we developed a new approach that address the batch effect issue in the preprocessing stage, resulting in better peak detection, alignment and quantification. It can be combined with down-stream batch effect correction methods to further correct for between-batch intensity differences. The method is implemented in the existing workflow of the apLCMS platform. Analyzing data with multiple batches, both generated from standardized quality control (QC) plasma samples and from real biological studies, the new method resulted in feature tables with better consistency, as well as better down-stream analysis results. The method can be a useful addition to the tools available for large studies involving multiple batches. The method is available as part of the apLCMS package. Download link and instructions are at https://mypage.cuhk.edu.cn/academics/yutianwei/apLCMS/

Metabolome Wide Association Study of serum DDT and DDE in Pregnancy and Early Postpartum

The advancement of high-resolution metabolomics (HRM) and metabolome-wide-association study (MWAS) enables the readout of environmental effects in human specimens. We used HRM to understand DDT-induced alterations of in utero environment and potential health effects. Endogenous metabolites were measured in 397 maternal perinatal serum samples collected during 1959-1967 in the Child Health and Development Studies (CHDS) and in 16 maternal postnatal serum samples of mice treated with or without DDT. MWAS was performed to assess associations between metabolites and p,p'-DDT, o,p'-DDT and p,p'-DDE levels, followed by pathway analysis. Distinct metabolic profiles were found with p,p'-DDT and p,p'-DDE. Amino acids such arginine had a strong association with p,p'-DDT and o,p'-DDT in both women and mice, whereas lipids and acyl-carnitine intermediates were found exclusively associated with p,p'-DDE in CHDS women indicating mitochondrial impairment. It suggests that the role of serine and fatty acid metabolism on the causal disease pathway should be examined

Exposome epidemiology for suspect environmental chemical exposures during pregnancy linked to subsequent breast cancer diagnosis

Breast cancer is now the most common cancer globally, accounting for 12% of all new annual cancer cases worldwide. Despite epidemiologic studies having established a number of risk factors, knowledge of chemical exposure risks is limited to a relatively small number of chemicals. In this exposome research study, we used non-targeted, high-resolution mass spectrometry of pregnancy cohort biospecimens in the Child Health and Development Studies to test for associations with breast cancer identified via the California Cancer Registry. Second and third trimester archival samples were analyzed from 182 women who subsequently developed breast cancer and 384 randomly selected women who did not develop breast cancer. Environmental chemicals were annotated with the Toxin and Toxin-Target Database for chemical signals that were higher in breast cancer cases and used with an exposome epidemiology analytic framework to identify suspect chemicals and associated metabolic networks. Network and pathway enrichment analyses showed consistent linkage in both second and third trimesters to inflammation pathways, including linoleate, arachidonic acid and prostaglandins, and identified new suspect environmental chemicals associated with breast cancer, i.e., an N-substituted piperidine insecticide and a common commercial product, 2,4-dinitrophenol, linked to variations in amino acid and nucleotide pathways in second trimester and benzo[a]carbazole and a benzoate derivative linked to glycan and amino sugar metabolism in third trimester. The results identify new suspect environmental chemical risk factors for breast cancer and provide an exposome epidemiology framework for discovery of suspect environmental chemicals and potential mechanistic associations with breast cancer